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Abstract 



It is shown how a recent method to systematically extrapolate and resum 
the loop expansion for nonlinear sigma-models is related to solutions of the 
renormalization group equation. This relation is used to generalize the explicit 
equations of state obtained previously to models which display crossover phe- 
nomena. As an example we discuss Wegner's localization model and consider 
the crossover from symplectic to unitary symmetry. 
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I. INTRODUCTION 



Nonlinear a-models have proven useful in Statistical Mechanics for the description of 
both ferromagnetism |1] and Anderson localization of noninteracting electrons The 
description by a nonlinear a-model allows to study these phase transitions in the vicinity of 
the lower critical dimension, which is d~ = 2 in either case. For the localization problem 
this description is particularly valuable, since no upper critical dimension is known for 
the localization, or Anderson, transition, no mean-field theory exists, and hence the usual 
methods of expanding about mean-field theory and the upper critical dimension are not 
applicable ||. The question then arises whether it is possible to derive a simple equation 
of state || for the localization problem, which is capable of qualitatively describing both 
phases and the transition between them. The answer turned out to be affirmative, but in 
a somewhat restricted sense. An approximate equation of state was derived for the case 
of localization in the presence of time reversal invariance |J. It turned out to correctly 
reproduce the critical behavior known from renormalization group (RG) methods to first 
order in an e-expansion about d = 2 [Q. This equation of state became known as the 
'self-consistent theory of Anderson localization', and it has been derived by a variety of 
techniques ||[7|]. However, it was not obvious how to generalize this theory to models with 
broken time reversal invariance, or how to systematically include higher orders in the e- 
expansion. Attempts to improve on this situation by exponentiating perturbation theory 
gave results that again were satisfactory only in the case of time reversal invariance ||. 

In a recent paper ||, to be denoted by I, the present authors have derived equations 
of state for nonlinear a-models that are based on a systematic loop expansion, and give 
physical results for localization models both with and without time reversal invariance. The 
method was also applied to nonlinear a-models for Heisenberg ferromagnets. The basic 
idea in I was to obtain the RG /5-function to a given order in the loop expansion, truncate 
it, and then construct a perturbation series for the relevant physical quantities (e.g. the 
conductivity in the case of localization, or the magnetization in the case of magnets) that 
yields the truncated /3-function exactly. The resummation of this extrapolated perturbation 
series led to transcendental equations whose solution gave the desired equation of state. 

Such an extrapolation and resummation of perturbation theory is essentially what the RG 
equation also achieves 0. Indeed, one of the derivations of the 'self-consistent localization 
theory' made use of the RG equation 0. It is therefore natural to ask whether there is a 
simple relation between these two methods. In the present paper we show that a simple 
relation does indeed exist, and how the RG equation method can be used to derive equations 
of state identical to those obtained in Ref . |§ . We then use this insight to derive an explicit 
equation of state for the case of localization with weakly broken time reversal invariance, 
when there is a crossover between two fixed points governed by different symmetries. The 
same method can be used to describe anisotropic Heisenberg ferromagnets [TIJ. However, 



the localization case is slightly more difficult since there the two fixed points governed by 
different symmetries appear at different orders in the loop expansion. Here we restrict 
ourselves to the localization problem, and just mention that essentially the same discussion, 
with some simplifications, applies to anisotropic magnets. 
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II. THE MODEL, AND ITS RENORMALIZATION 



A. Definition of the Model 

As in I, we consider Wegner's model for the localization of noninteracting electrons 0. 
Both the model and our notation will be the same as in I, and so we will be very brief. 
Wegner's model is a nonlinear cr-model defined on the symmetric spaces Sp(2N) / Sp(N) x 
Sp(N) and U(2N)/U(N) x U(N) for systems with and without time reversal invariance, 
respectively. Quenched disorder is handled in these models by considering N replicas of the 
system, and taking the limit iV — > ("replica limit") after calculations. The action reads 

S = ^Jdxtr (VQ(x)) 2 + ^Jdxtr [AQ(x)] , (2.1a) 

where Q is a 2N x 2N matrix subject to the constraints, 

Q 2 (x) = 1 2JV , tr Q(x) = , (2.1b) 



and A is a block-diagonal matrix, 



-J ■ 

with 11 jv the N x N unit matrix. The symmetry of the action is symplectic or unitary, 
respectively, depending on whether the matrix elements of Q are quaternion or complex 
valued. We will refer to these models as the symplectic and the unitary model, respectively. 
In both cases the matrix Q is hermitian, and in the symplectic model it is subject to the 
additional constraint Q = C T Q T C, where C — ia x ® a y with Pauli matrices a XjV . 

It is well known that each of the two models defined by Eqs. (|2.1|) describes a localization 
transition of noninteracting electrons fj,ll|- The bare coupling constant T in Eq. (|2.1a ) 



is proportional to the conductivity in self-consistent Born approximation (SCBA); it is a 
measure of disorder. H = QnNp/A with Q the external frequency, and Np the density 
of states at the Fermi level. The symplectic model describes the localization transition 
in systems that are time reversal invariant. If time reversal invariance is broken, e.g. by 
magnetic impurities, the critical behavior is instead described by the unitary model. 



For later reference we also define an additional term [12fl , which we add to the action for 
the symplectic model, 

G 



6S= 4T J d * tr [Q(X) ' ~ iUz ® l7v]2 ' (2 ' 2) 



where [a, b] denotes the commutator of a and b. The term SS breaks the symmetry of the 
symplectic model action, and for any nonzero G the asymptotic critical behavior of the 
phase transition described by the symplectic model plus the term 5S is the same as that 
of the unitary model. For small values of G, however, the asymptotic critical behavior is 
preceded by a region that displays a crossover between the symplectic and the unitary fixed 
point. This model is relevant for localization with weakly broken time reversal invariance, 



e.g. in the presence of a small concentration of magnetic impurities [12]. It also describes 
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the localization of third sound waves in superfTuid He films on a rough substrate with a 
superimposed uniform superfTuid flow [|i"3H . 

It is convenient to eliminate the nonlinear constraint, Eq. ( |2.1b| ), by parameterizing the 
matrix Q as, 



Q 



H- N - qq 



1/2 



N 



Q 



t 



q 

ILjv - q ] q 



1/2 



+ a 



(2.3) 



N , 



and to expand q and q^ in a quaternion basis, 



q = J2 1r 



(2.4) 



where r = 0,1,2,3 and r = 0, 3 for the symplectic and unitary models, respectively, and 
t = U 2 , t 1j2) 3 = —icr x ,y,z- The g r are N x N matrices. 

Let us consider the two-point q- vertex function, which is the inverse of the q-q propagator. 
The most general model is the symplectic one with the symmetry breaking term, Eq. ( |2.2|) , 
added. One obtains two different 2-point vertex functions for the channels r = 0, 3 and 
r = 1,2, respectively. To Gaussian order they read, 



r«2*(p) 



i 

T 



V 



H 



O T 



(2.5a) 



T^ p (p) = ^(p 2 + H + G)+0(T°) . (2.5b) 

The subscripts p — h and p — p refer to the particle-hole (r = 0,3) and particle-particle 
(r = 1,2) channels, respectively |j4,TT|] . G = corresponds to the symplectic model. In 



this case both vertex functions coincide and are soft, or massless, i.e. they vanish with 
vanishing wavenumber and frequency. The soft modes describe diffusion in the particle- 
hole and particle-particle channels. A nonzero G breaks time reversal invariance, and gives 
the particle-particle channel a mass. The particle-hole channel remains massless because of 
particle number conservation. 

The most interesting physical quantity in the context of the localization problem is the 
conductivity, a. It is related to T^_ h by 

a = H& Tr ™ M ' (2 ' 6) 



where we have normalized a by its SCBA value. 



B. Renormalization 

It is well known that the symplectic and unitary models defined by Eqs. (|2.1| ) are renor- 
malizable with two renormalization constants, one for the coupling constant T, and one field 
renormalization constant 0] . In the replica limit the field renormalization vanishes, and one 
has only one nontrivial renormalization constant, which we denote by Z t . The symmetry 
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breaking term, Eq. (|2.2|) , requires one additional renormalization constant, Z g) for the cou- 
pling constant G. Denoting the renormalized disorder, frequency, mass, and conductivity by 
t, u, g, and ay, respectively, the relations between the bare quantities and their renormalized 
counterparts are 



T = K~ £ Z t t , 


(2.7a) 


Q = K 2 Z t u , 


(2.7b) 


G = K 2 Z g Z t g , 


(2.7c) 


a = Z t a r , 


(2.7d) 



with e = d — 2, and k an arbitrary RG momentum scale. From the requirement that the 
bare theory must be independent of k one obtains a RG equation || for oy, 



- t P(t, g) + P{t, g) | + y(t,g)§- g + (i«t, S ) - d) U A 



<7 r (t,w, 9 ) = . (2.8) 



Here /3(t,g) = dt/d Inn and j(t,g) = dt/d\nK, with the derivatives taken at a fixed bare 
theory, are RG functions. They can be obtained in a loop expansion from the perturbation 
expansion for the vertex function, Eqs. (|2.5| ), by standard methods. We have performed the 
calculation to two- loop order with the result, 

(3(t, g) = et - - — — — 9 — + 0(t 4 ) , (2.9a) 

MV ,yj 4 1 + # 32 1 + # V 7 ' V ; 

1 (t,g) = -2g + - -^- + 0{e) . (2.9b) 
4 1 + 3 

In the limits g = (symplectic model), and g = oo (unitary model) Eq. ( ^.9a| ) reproduces 
well-known results [BJ. The /3-functions for these two cases are actually known to five-loop 
order |14 , [15| . However, for the sake of simplicity we will restrict ourselves to the two- loop 



result as contained in Eqs. ( |2.9| ). Consequently, all critical exponents, scaling functions, etc. 
will be valid to two-loop order only. Inspection of Eqs. (|2.9| ) shows that they allow for two 
nontrivial fixed points (FP), viz. (t*,g*) = (4e,0), and (t*,g*) = (y/32e,oo). The former 
is unstable with respect to g, while the latter is stable. These two FP correspond to the 
symplectic and unitary universality classes, respectively, of the Anderson transition. 



III. THE EQUATION OF STATE 

We now obtain the desired equation of state, i.e. the conductivity as a function of 
disorder t, frequency u, and symmetry breaking parameter g, from the RG equation, Eq. 
( p.8[ ). For simplicity we do this for the renormalized theory. Since all bare and renormalized 
quantities are proportional to each other, the same functional relations hold for both the 
bare and the renormalized theory, and we will drop the subscript on oy. 
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A. The Symplectic and Unitary Cases Revisited 



Let us first reconsider the symplectic and unitary models. In either case there is no 
(^-dependence, and the /^-function is given by Eq. ( |2.9a| ) with g = or g = oo, respec- 
tively. At zero frequency the PDE ( |2.8|) (without the d/dg-tevm) then turns into an ODE, 
which is readily solved in terms of a quadrature. With our normalization condition for the 
conductivity, a(t — > 0,u = 0) — > 1 (see Eq. (|2.6|) , we obtain for the static conductivity, 
a (t) = a(t,u = 0), 



<Jo(t) = exp / dr 



(3(i 



(3.1) 



Of course, cxo(t) = is also a solution of the RG equation. The physical conductivity 
is determined by Eq. ( |3.1|) only for t < t c , where t c is given by j3{t c ) = 0. For t > t c the 
physical conductivity vanishes, and o"o(t) as given by Eq. ( |3.1|) , and interpreted as a principal 
value integral, has no direct physical meaning. For t < t c Eq. ( |3.1| ) yields, 



a (t) =c(t)\l-t/t c 



(3.2a) 



where the function c(t) is analytic at t — t c , and the conductivity exponent s is related to 
the correlation length exponent v, 



ve 



-e/(3'(t c 

The first equality in Eq. ( |3.2b[ ) is known as Wegner's scaling law [ IE ] 



(3.2b) 



Using standard methods JT7j it is now easy to express the general solution of the PDE, 
Eq. ( |2.8|) , for g = or g = oo in terms of a {t). The associated system of ODEs is 



d t 



d In a; 



d In a 



P(t) -d + /3(t)/t -e + /3(t)/t 



(3.3) 



Let us consider the equalities relating the first term and the third term, and the first term 
and the second term, respectively. The general solutions of these two ODEs are, 



«i = cr/a (t) 



(3.4a) 



Bj^^Wt)]-*' , (3.4b) 

where u\ and U2 are integration constants. The general solution of the PDE is F(u\ 1 U2) = 0, 
with F an arbitrary function of two variables. This can be rewritten as, 

a(t,u) = a (t) F(ut 2 ^ao d/e (t)) , (3.5) 

with F an arbitrary function of one variable. In order to determine F we need to impose 
a boundary condition. Since the conductivity at large frequencies (u ~ 1) is qualitatively 
insensitive to the disorder, we can normalize the conductivity by its value at uj = 1, 

a(t,cu = l) = l . (3.6) 
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From Eq. ( |3.5|) this yields a transcendental equation for F, which can be rewritten as 



x = [t F (F(x))} 2 ^ [F(x)f* , 
where tp(x) is obtained as the inverse of the relation between F and a , 

F=l/a {t) -> t = t F {F) . 



(3.7) 



(3.8) 



Note that tp(F) has two branches. Finally, from Eq. ( |3.5| ) we know that F(x) = cr/cr (t) 
with the argument of the function F given by x = tut 2 ^ [<7 (t)] _d//<E . Using this is Eq. 
we obtain the equation of state in the form, 



Zdfe 



(3.9) 



Equation ( |3.9| ) is identical with the result that was obtained by different methods in I. 
This shows that the resummation scheme of I is equivalent to solving the RG equation, Eq. 
(|2.8|) , with the boundary condition given by Eq. (|3.6| ). Let us write down explicit results by 
using the /5-function, Eq. ( |2.9a| ), to two-loop order. The generalization to higher loop orders 
is straightforward, and has been given in I. To two-loop order we obtain for the symplectic 
model, 



-2/e 



with t r 



a d/e 

4e, and for the unitary one, 



t/t c 



a 



2/6 



(3.10a) 



CO 



a 



d/e 



(t/t c 



,-2/e 



(3.10b) 



with t c = V32e. Equation ( 3.10a ) is the well-known result of Ref. 0. The above derivation is 
essentially identical with the one given by Hikami 0. Equation (|3.10b| ) is the corresponding 
result for the unitary model, which was first derived in I. We note that the derivation using 
the RG equation is shorter and more transparent than the resummation scheme used in I. 



B. The Crossover from Symplectic to Unitary Critical Behavior 

We will now apply the same method to construct an equation of state for our more 
general model, defined by the symplectic model with the symmetry breaking term, Eq. 
The RG equation is given by Eqs. (|2.8|), fl2.9p. Let us first consider the zero frequency 



conductivity, a (t,g) = a(t,uj = 0,g). In this case the system of associated ODEs consists 
of two independent equations, which we choose as, 

dt d9 (3.11a) 



P(t,g) j(t,g) ' 

d hio-p d g miM 

(3(t,g)/t-e j(t,g) l ' ' 
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The general solution of Eq. ( |3.11a|) can be written in the form, 



f 1 {t,g)=ui , 

with U\ an integration constant. The general solution of Eq. ( 3.11b ) reads, 

/3(t,g)/t' 



ctq exp 



dg- 



«2 



(3.12a) 



(3.12b) 



where U2 is a second integration constant. The general solution of the PDE for o"o is given 
by F(ui,U2) = with F an arbitrary function, or 



<ro(t,g) = /(ui) a (t,g) , 



(3.12c) 



with Ui a function of t and g according to Eq. ( p,12a| ), / an arbitrary function of Ui, and 

P(t,g)/t-e\ 



&o(t,g) = exp f / dg- 



i(t,g) 



(3.12d) 



Note that for determining the indefinite integral in Eqs. ( |3.12b| ), ( |3.12d[ ), t must be expressed 
in terms of g and u± by means of Eq. ( |3.12a|) . U\ is kept fixed for the integration, and then 
replaced by fi(t,g) again. 

Equation fl3.12c|) is the desired equation of state at zero frequency. What remains to 
be done is to determine u% or f\ explicitly as a function of t and g, and to determine the 
unknown function / by imposing suitable boundary conditions. This is a substantially more 
complicated task than for the pure symplectic or unitary model. We first derive an explicit 
expression for the function fx, Eq. ( |3.12a| ). Even with the simple two-loop approximations 
for the RG functions, Eqs. (|2.9|) , we have been unable to find a closed form representation 
for this function. However, it can be discussed analytically in certain limits. Let us first 
consider the case of weak symmetry breaking, g « 1. For g = there is a fixed point 
already at one-loop order, and we can restrict ourselves to a one-loop approximation for the 
/3-function, j3(t,g) « et — (t 2 /4)(l — g). Since both g and t are small compared to unity, it 
is then sufficient to use the zero-loop approximation for the RG function 7, "f(t,g) ~ —2g. 
Equation ( |3.11a ) is then readily integrated with the result, 



h(t,g^o) 



g 



e/2 



(3.13a) 



In the same limit one obtains from Eq. ( |3.12d| ) <7o(i, g) = tg t l 2 . With t° c = t c (g = 0) = 4e the 
critical disorder for the symplectic model, the zero frequency conductivity can be written, 



a (t,g^0)=tg^ 2 f 



, tg^ , 



(3.13b) 



For t near t® we thus obtain the scaling form which was first derived by Khmel'nitskii and 
Larkin from a scaling hypothesis We see that both the value of the conductivity at 

t = t® and the shift of the mobility edge are proportional to g 1 ^ with <\> = 2/e + 0(1) the 
crossover exponent |L8| , |l9|l . Comparing Eqs. ( |3.13b| ) and ( |3.10a|) we also see that for large 
arguments (i.e., for g — > 0), the function / behaves like f(x — > 00) ~ x. 
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Similarly, one obtains in the limit of large g, 



h{t,g-> oo) = — 



1 1 1 

+ 



1 


"1 1 


= ? 


i* 2 t*(g)\ 



t 2 32e 32g 

For ao one obtains again cr (i, g) = £<p 2 , so we have for the conductivity, 



(3.14a) 



(3.14b) 



Close to criticality this leads to the scaling form 



1/0 



(3.14c) 



with 0' = — 1/e. Here we have used that /(a; — > 0) ~ y^F, which can be seen by comparing 
Eqs. ( p. 14 b| ) and ( p. 10b ). Notice the qualitative difference between the scaling forms at small 
g, Eq. ( p,13b|) , and at large g, Eq. (|3.14c| ). In the former case, the conductivity is driven 
to zero because of a zero in the scaling function /, and the asymptotic critical behavior is 
determined by the behavior of / in the vicinity of that zero. This gives rise to a crossover 
from critical behavior governed by the symplectic fixed point to the asymptotic behavior 
which is governed by the unitary fixed point. We will see this crossover behavior in more 
detail below, when we determine the function / explicitly. In the latter case the conductivity 
vanishes due to the prefactor in Eq. ( 13.14a) ). The function $ merely provides corrections to 
this asymptotic scaling behavior, and there is no crossover. The exponent 0' is a corrections 
to scaling exponent. The reason for this difference is that g is a relevant operator with 
respect to the symplectic fixed point, while 1/g is an irrelevant operator with respect to the 
unitary one. It is well known that a relevant operator leads to a crossover, while irrelevant 
ones lead to corrections to scaling. Our RG treatment thus does not confirm a previous 
conjecture that the scaling form of the conductivity is the same for both small and large g 
0- _ 

We now turn to the task of determining the scaling function fi for all g. A problem one 
encounters in that context is that the solution of the equation fi(t, g) = 0, which determines 
the critical disorder of the system as a function of g, develops from a single root, t = t° c for 
(7 = 0, into a double root, t = ±t£° for g = oo. This problem can be solved as follows. It 
is convenient to switch from the variables (t,g) to new variables (x, y) = (l/t,g/(l + g)). 



Equation (|3.11a|) reads in terms of these variables, 

l-y)/A + y/32x 



dx 
dy 



-ex 



-2y{l-y)+y\l-y)/Ax 



(3.15) 



Let us now consider flow equations for x and y that lead to the same fx, determined by Eq. 
( |3.15| ), as the original ones, but are obtained by multiplying numerator and denominator on 
the r.h.s. of Eq. ( |3.15|) , by —x, 



dx 
~dl 



ex 



ar(l-y)/4-y/32 



(3.16a) 
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4 = 2xy(l -y)-y 2 {l-y)/A 



(3.16b) 



These flow equations allow for four fixed points, viz. (#1,2/1) = (l/4e, 0), (22,2/2) = 
(l/v32e, 1), (#3,2/3) = (— 1/V32e, 1), and (24,2/4) = (0,0), and the flow diagram is shown 
in Fig. 1. Note that only the region x > 0,2/ > in this diagram is physical, and that 
the region x < serves a purely mathematical purpose. Important trajectories in this flow 
diagram are the curves x = x+{y)i an d x = X-{y)i which connect the fixed points (1) and 
(2), and (3) and (4), respectively. We have determined only the qualitative behavior of these 
curves to the extent that we will need them, viz. 

X+(y -»• 0) = #1 - const x 2/ V * , (3.17a) 



X +(y^l) = x 2 + Ul + 0(Ve))(l-y) , 



and 



X-(2/-0) = 



(3.17b) 



(3.17c) 



x _(y_>l) = _s a + i(l + 0(Vi))(l-y) . 



(3.17d) 



The constant in Eq. [3.17a| can only be determined numerically, and of X-{u 0) we wm 
need to know only that it vanishes. We now write fx as a function of x and y in the following 
form, 



y- e / 2 (\ - yY 

/i(#,2/) = [x-x+(y)\ [x-x-(y)\ y + x ( 1 + y y/2 G ( x >y) 



(3.18) 



The explicit part of this equation contains all qualitative features of the flow. In particular, 
all of the singular behavior has been built in, as well as the known behavior at small and 
large g, respectively. The function G(x,y) is therefore a harmless function of both of its 
variables, and we approximate G(x,y) ~ 1. Changing the variables back to t and g, we then 
have 



f 1 (t,g) = [l/t-l/t*(g)} [l/t-l/t c {g)\ 



9 1 

+ 



1 



e/2' 



1 + 2? t Vl + o 



where 



t c (g — > 0) = 4e + const x g 



e/2 



g-^(l+g)-^ , 

(3.19a) 
(3.19b) 



tc(9 00) = v32e — const/5 1 , 
is the critical disorder as a function of g, and 



(3.19c) 
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t*(g ^ 0) ^ oo , (3.19d) 

t*{g -»• oo) = -V32^- const/# . (3.19e) 

is the auxiliary curve corresponding to X-(y)- 

Now we are in a position to construct an approximate equation of state for the case of 
zero frequency that contains all qualitative features which we have discussed so far. Let us 
return to Eq. Q3.12c|) . As we have seen, 5"o is equal to tg e ^ 2 both for small g and for large 
g, and we approximately put <To = tg e ^ 2 everywhere. From Eqs. ( 3.13b ) and (|3.14b| ) we see 



that the function f(ui) is proportional to U\ for large U\ and proportional to ^/ui for small 



U\. This function can therefore be adequately modeled by f(x) = yx(l + x). We now have, 

a Q {t,g)=tg^ 2 \fi(t,g)(l + hfag))] 1 ' 2 , (3.20) 



with fi(t,g) from Eqs. ( gT§ 



Equations ( |3.19|) , ( |3.20|) are the desired result for the equation of state, i.e. for the 
conductivity as a function of the disorder t and the symmetry breaking parameter g at zero 
frequency. Let us briefly discuss this result. For g — » 0, 1/t* vanishes, t c approaches 4e, 
fi diverges like g~ e ^ 2 , and a becomes independent of g. Equation ( |3.20| ) then reproduces 
the result for the symplectic model, Eq. (|3.10a|) at u = 0. For g — > oo, 1/t* approaches 
—l/t c (g = oo) = — 1/V32e, and fi vanishes like g~ e . Again cr becomes independent of g, 
and Eq. (|3.20|) reproduces the result for the unitary model, Eq. ( |3.10b|) at uj = 0. If one 
approaches the critical disorder, t — > t c (g), at fixed but small g (g « 1), then the critical 
region consists of two regimes. For |t — t c (g)\ >> g e ^ 2 one has fi >> 1, and <To shows a 
behavior characteristic of approaching the symplectic localization transition. In particular, 
the critical exponent s = 1. However, when \t — t c (g)\ « g e ^ 2 one has /i << 1, and 
the asymptotic approach to the phase transition is governed by the critical behavior of the 
unitary model. This is the crossover behavior induced by the instability of the symplectic 
fixed point with respect to the relevant operator g. With increasing g the crossover region 
moves away from the critical point, and for sufficiently large values of g it is outside of 
the critical region. The critical behavior is then entirely that of the unitary model, with 
corrections to scaling due to the irrelevant variable 1/g. 

We finally take the frequency into account. We choose the third ODE associated with 
the PDE, Eq. (gj), as, 

dlnU d9 . (3.21a) 



-d + p(t,g)/t -y(t,g) 
The general solution reads, 

ut 2 ^[a (t,g)]- d ^ = u 3 , (3.21b) 

with u 3 a third integration constant, in addition to ui and u<i- The general solution of the 
full PDE can then be written as, 

a(t, oo, g) = a (t, g) G (f^t, g),ujt 2 ^ [a (t, g)}~ d/e ) , (3.22) 
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with G an arbitrary function of two variables. The function / of Eq. (|3.12cj) is related to G 
by G{x, 0) = f(x). We now impose the same boundary consition as in Sec. [Ill A , namely 



a(t,u = l,g) = 1. The transcendental equation for G(ux, x) analogous to Eq. (|3.7| ) then 
reads, 

x = [t G (G{u l7 x), Ul )} 2/t [G( Ul , x)} d/e , (3.23) 

Here t G , which is a function of two variables, is obtained as follows. In the relation G = 
l/ao(t,g) we express g = g(t,u{) as a function of t and u\ according to Eq. ( |3.12a| ). t G is 
then given explicitly function of G and u\, 

G = l/a {t,g{t, Ul )) -> t = t G (G, Ul ) . (3.24) 

Now we restore the frequency dependence. With x = ut 2 / e [ao(t, g)\~ d ^ one has G(ui,x) = 
a /gq. Inserted into Eq. ( ^.23|) this yields, 

X t- 2 IU G {a/a Q (t,g),h(t,g)) . (3.25) 



a d/e 

This is the desired equation of state for the general model. It contains the previous results 
for the symplectic (g = 0) and unitary (g = oo) models as special cases, and explicitly 
describes the crossover as demonstrated above. It also gives a physically correct description 
of both the metallic and the insulating phases. This has been discussed in detail in I, and 
can be seen from Eqs. ( |3.1U| ), in the limits g = and g = oo. Since the frequency dependence 
of the conducticity is qualitatively the same in both limits, this is also true for arbitrary 
values of g, and hence for Eq. ( |3.25| ). 
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FIGURES 

FIG. 1. Schematic flow diagram for the flow equations, Eqs. ( |3.16| ) 
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